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We have studied solid hydrogen up to pressures of 300 GPa and temperatures of 350 K using 
density functional theory methods and have found "mixed structures" that are more stable than 
those predicted earlier. Mixed structures consist of alternate layers of strongly bonded molecules 
and weakly bonded graphene-like sheets. Quasi-harmonic vibrational calculations show that mixed 
structures are the most stable at room temperature over the pressure range 250-295 GPa. These 
structures are stabilized with respect to strongly-bonded molecular phases at room temperature by 
the presence of lower- frequency vibrational modes arising from the graphene-like sheets. Our results 
for the mixed structures are consistent with the experimental Raman data [M. I. Eremets and I. A. 
Troyan, Nature Mater. 10 927 (2011) and R. T. Howie et al. Phys. Rev. Lett. 108 125501 (2012)]. 
We find that mixed phases are reasonable structural models for phase IV of hydrogen. 



I. INTRODUCTION 

Hydrogen is the most abundant and ancient element in 
the Universe and the simplest of all atoms. Hydrogen was 
solidified by Dewar— , and Wigner and Huntington pre- 
dicted that solid hydrogen would metallize at pressures 
of about 25 GPa£, although the metallization pressure 
is now known to be much higher. Insulating molecular 
phases of solid hydrogen have been identified in static 
diamond anvil cell (DAC) experiments up to pressures 
of over 300 GPa. Phase I is a quantum crystal con- 
sisting of rotating hydrogen molecules on a hexagonal 
close packed lattice^. The orientationally ordered molec- 
ular phases II and III^ appear at pressures of about 110 
and 150 GPa, respectively^, but their structures have 
not been determined experimentally Ashcroft suggested 
that metallic hydrogen might superconduct up to room 
temperature^ and, more recently, possible superconduct- 
ing superfluids and metallic superfluids of hydrogen have 
been suggested^. The gas giants Jupiter and Saturn and 
some of the recently discovered extrasolar planets^ are 
believed to contain large amounts of metallic hydrogen 
at high temperatures and pressures into the terapascal 
range. Metallic hydrogen has been generated at temper- 
atures of several thousand Kelvin in dynamic shock wave 
experiments 8 -, but it has proved very difficult to form 
metallic hydrogen at lower temperatures in DACs. 

DAC studies at pressures up to 300 GPa are approach- 
ing current experimental limits. Determining crystalline 
structures of hydrogen at high pressures is extremely dif- 
ficult because x-ray diffraction, which is the workhorse of 
crystal structure determination, is severely hampered by 



the weak scattering of hydroger 



suo 



Rama: 



.3.11.12 



and 



infra-red (IR^ii 3 . vibrational spectroscopies have been 
very important in studying the bonding and structures 
of hydrogen at high pressures, although they have not 
proved sufficient to determine the structures themselves. 

Recent room temperature (295-300 K) DAC experi- 
ments up to pressures of around 300 GPa have found 



evidence for a previously unreported structural phase 
of hydrogen, that we refer to as the putative phase 
I V 14 i 15 . The intense Raman molecular vibron mode 
started to drop rapidly in frequency as the pressure was 
increased above about 200 GPa, and at higher pressures 
the samples became dar k 14 i 15 . Darkening of samples 
above 200 GPa was observed in earlier low-temperature 
experiments 1 ^ 1 ^, suggesting a reduction in the band gap 
that may be a precursor to metallization. Darke ning 
was also observed at high pressures in Refs. GJ and (Tf 
Eremets and Troyan reported an increase in electrical 
conductance around 260-270 GPaM. The Raman data 
in Refs. [TJ andll5l are similar, although in Ref. [TH signals 
characteristic of phase IV were obtained up to higher 
pressures of 315 GPa. The Raman data of Refs. [13 
and [TH can be used to eliminate candidate structures for 
phase IV. There arc differences between these experimen- 
tal data, but here we focus on quantities for which there 
is a reasonable degree of consistency between the studies. 

In previous wor k 17 ' 18 we searched for high-pressure 
structures of hydrogen using density functional theory 
(DFT) methods and ab initio random structure search- 
ing (AIRSS) 1 ^^. Structures were found at pressures 
from about 80 to nearly 400 GPa which were calculated 
to be more stable than previously-known ones. Several 
candidates for phase II were found, although the en- 
thalpy differences between them are very small. A single 
outstanding candidate of C2/c symmetry was found for 
phase Illi 7 -. This structure was the lowest in enthalpy 
from about 70-250 GPa when proton zero-point (ZP) 
motion was included, which covers most of the region in 
which phase III is observed. The C2/c structure also ex- 
hibits strong IR activity that is characteristic of phase 
Il£L2>2i, and it gives a good account of the Raman-active 
vibron frequency. A molecular structure of Cmca sym- 
metry was found to be stable in the range 250-380 GPa, 
which we refer to as Cmca-12, where the "12" denotes 
the number of atoms in the primitive unit cell. Above 
380 GPa another structure of Cmca symmetry^ was the 
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most stable, which we refer to as Cmca-A. The discov- 
ery of a class of "mixed" hydrogen structures consist- 
ing of alternate layers of strongly bonded molecules and 
weakly-bonded graphene-like sheets was reported in Rcf. 
[l7l . Mixed structures with Pbcn, C2 and Ibam sym- 
metries were found and, when ZP motion was included, 
Pbcn in particular became energetically competitive with 
the C2/c phase over a wide range of pressures 1 ^. 

In this work we have used AIRSS to identify candi- 
date structures for phase IV. In the AIRSS approach 
randomly chosen structures are relaxed to minima in the 
enthalp y 19 ' 20 . The basic AIRSS method has virtually 
no free parameters and is essentially unbiased, and it is 
therefore an ideal starting point from which to introduce 
useful biases and constraints. The most important tech- 
nical improvements on the searches of Ref. [l?] and [TH 
are that here we employ constraints to exclude very-low- 
symmetry structures which, according to Pauling's rule 
of parsimony 2 ^, are unlikely to be low in enthalpy, and 
that we have "shaken" low-enthalpy structures by mak- 
ing small random atomic displacements and re-relaxing. 
We have also shaken supercells of low-enthalpy struc- 
tures that are then able to relax into structures with 
the larger supercell periodicity. AIRSS has been ap- 
plied to many systems including molecular solids such as 
hydroge n 17 ' 24 , water—, nitrogen 2 ^, ammonia^ 7 - and am- 
monia monohydrate 2 *. 



II. DESCRIPTION OF THE CALCULATIONS 
A. DFT calculations 

Our main DFT calculations were performed with the 
CASTEP 2 ^ code, the Perdew-Burke-Ernzerhof (PBE)^ - 
Generalized Gradient Approximation (GGA) density 
functional, and ultrasoft pseudopotentials^i. For the 
searches we used a standard ultrasoft pseudopotential, 
a plane wave cut-off energy of 230 eV and a Monkhorst- 
Pack^ Brillouin zone sampling grid of spacing 2ir x 0.05 
A -1 . The structures obtained in the searches were re- 
relaxed at a higher level of accuracy consisting of a hard 
ultrasoft pseudopotential, a plane wave cut-off energy of 
1000 eV and a Brillouin zone sampling grid of spacing 
2tt x 0.03 A" 1 . The enthalpies of the structures reported 
in the paper were calculated at the higher level of ac- 
curacy. For the phonon, Raman and IR calculations we 
used a plane wave cut-off energy of 1200 eV and Bril- 
louin zone sampling grids of spacing 2n x 0.04 A^ 1 . The 
phonon calculations were performed with the CASTEP 
code and cross-checked with the Quantum ESPRESSO 
code2£. 



B. AIRSS calculations 

Wc have performed AIRSS runs in addition to those re- 
ported in our earlier paper s 17 ' 18 . The searches were per- 



formed at 250 GPa unless otherwise specified. Searches 
in 24 and 48 atom cells using symmetry constraints with 
from 2 up to 12 symmetry operations were performed, 
and 4956 relaxed structures were obtained. These struc- 
tures included the important Pc structure (Fig. [I]) which 
was found in searches with 2 symmetry operations. The 
16-atom Ibam mixed structure^ 7 - was shaken in super- 
cells with up to 96 atoms, and 305 relaxed structures 
were obtained. The 48-atom Pbcn mixed structure was 
shaken in supercells with up to 192 atoms, giving 156 re- 
laxed structures, and the Pc structure was also found in 
these searches. The 24-atom C2/c structure was shaken 
in supercells with up to 96 atoms and 253 relaxed struc- 
tures were obtained, and the 48-atom Pc structure was 
shaken in supercells with up to 192 atoms at 300 and 350 
GPa, giving 333 relaxed structures. We also found the 
Cmca- 12 and Cmca-A structures in our searches. Note 
that the C2/c structure that is our best candidate for 
phase III has 24 atoms in the primitive unit celU 7 - and 
that it is different from the C2 jc structure found by Tse 
et al£^, which has 12 atoms in the primitive unit cell. In 
the supplementary information for Ref. [13 our 24-atom 
C2/c structure was referred to by an alternative setting 
of B2/n. The two C2/c structures differ in the stacking 
of layers, and the 24-atom structure is lower in enthalpy. 
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FIG. 1. (color online). Side and top views of the Pc structure 
at 250 GPa. The side view shows the four-layer repeat with 
a weakly-bonded graphene-like layer on top (shaded atoms). 
The red lines show the next-nearest contacts within the layer. 



III. RESULTS 

A. Structure searches 

The searches confirmed our previous results 1 ^, but we 
also found the Pc symmetry mixed phase which has 48 
atoms in the primitive unit cell^. The Pc structure is 
similar to Pbcn, but with a somewhat larger distortion 
of the graphene-like layers, see Fig.[T] Static lattice DFT 
calculations show Pc to be about 2 meV per atom more 
stable than Pbcn, and at 200 GPa it is only about 3 meV 
per atom higher in enthalpy than C2/c. We emphasize 
that the enthalpy differences between competing phases 
are particularly small in hydrogen, with shifts of order 1 
meV per atom leading to significant changes in the nu- 
merical results. Howie et alJ£ recognized that the Raman 
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data for phase IV might be explained by a mixed struc- 
ture such as Pbcn, based on the calculated data reported 
in Ref . [3 



B. Bond lengths of the structures 

As shown in Fig. [2j the shortest bonds in C2/c and 
Cmca-Yl at 150 GPa of lengths 0.745 and 0.753 A, re- 
spectively, expand slowly with pressure, reaching 0.753 
and 0.758 A at 300 GPa. Fig.[T]shows that the molecules 
in the weakly bonded layers of Pc (and similarly for Pbcn 
and C2) form distorted hexagons. The structure evolves 
into the higher symmetry Ibam space group above about 
350 GPa, in which the weakly bonded layers adopt the 
graphene structure. The shortest bond in Pc at 150 GPa 
of length 0.729 A slowly contracts with increasing pres- 
sure, reaching 0.722 A at 300 GPa. The bond lengths in 
the weakly bonded layers of Pc, however, change more 
rapidly with pressure. At 150 GPa the weakly bonded 
layers have first and second nearest-neighbor distances of 
0.763 and 1.175 A, respectively, but upon further com- 
pression the shorter bonds expand and the longer ones 
contract, and at 300 GPa their bond lengths are 0.808 
and 0.962 A, respectively. These bond lengths become 
more nearly equal as the layer approaches the graphene 
structure. 
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C. Phonon calculations 

We performed quasi-harmonic vibrational calculations 
for the C2/c, Cmca-12, Cmca-4, Pbcn and Pc phases, 
see Fig. [3] The C2/c, Cmca-12 and Cmca-4 phases are 
dynamically stable at 250 GPa, while the Pbcn and Pc 
phases have some unstable modes. The unstable modes 
of Pbcn cover a significant portion of reciprocal space, 
while Pc is lower in enthalpy than Pbcn and has unstable 
modes over only a small region of reciprocal space. To 
investigate this further we doubled the unit cell of Pc to 
96 atoms and rc-relaxcd, obtaining a Pc-96 structure^ 
of the same symmetry that is slightly lower in enthalpy 
than Pc and has stable phonon modes (Fig. [3]). The 
large size and low symmetry of the Pc-96 structure make 
it difficult to calculate accurate vibrational spectra, and 
therefore we used the Pc phase as our model for phase 
IV. 



D. Enthalpy- pressure relations and phase diagram 

We calculated the proton ZP energy of the C2/c, 
Cmca-12, Cmca-4 and Pc structures at several volumes 
and constructed the enthalpy at zero temperature both 
with and without ZP motion, and the Gibbs free energy 
at 300 K, as shown in Fig.|3J At the static lattice level, Pc 
is less stable than C2/c over the pressure range covered 
by (Fig. 2k) and it becomes less competitive compared 
with C2/c and Cmca-12 as pressure is increased. The 
Cmca-A phase is more than 10 meV per proton higher in 
enthalpy than C2 Jc at these pressures, and it is therefore 
not visible in Fig. S^i). Including the vibrational ZP mo- 
tion (Fig. changes the picture significantly as the Pc 
and Cmca-4 structures become more energetically com- 
petitive because they have lower vibrational frequencies 
than the strongly molecular C2/c and Cmca-12 phases. 
Adding vibrational effects at a temperature of 300 K (Fig. 
Hfc) results in the Pc phase becoming the most stable in 
the range 245-300 GPa. The corresponding phase di- 
agram (Fig. [5]) shows a region of stability for Pc that 
is compatible with the experimental data for phase IV. 
Below 260 K, we also find a transition from the C2/c 
to Cmca-12 phases at about 285 GPa, which suggests 
the stability of another phase of molecular hydrogen at 
pressures beyond phase III. 



FIG. 2. (color online) Evolution of the bond lengths of C2/c, 
Cmca-12 and Pc with pressure. The molecular bond lengths 
of C2/c and Cmca-12 slowly expand with increasing pressure, 
while the short bonds in the strong layers of Pc undergo a 
slight contraction. The most notable feature is the behavior 
of the bonds in the weak layers. The short bonds in the weak 
layers expand with pressure rather more rapidly than those of 
C2/c and Cmca-12, while the long bonds in the weak layers 
contract rapidly with pressure. 



E. Raman spectra 

To investigate the possible relationship between phase 
IV and the most favorable structures obtained in our 
searches, we have calculated the frequencies and intensi- 
ties of the Raman modes of C2/c, Cmca-12, Cmca-4 and 
As can be seen in Figs. [5] and [7] and Table HI the 
C2/c and Cmca-12 phases each have a single intense Ra- 
man vibron mode whose frequency slowly declines with 
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FIG. 3. (color online) Phonon dispersion relations of (a) Pc, (b) Pc-96, (c) Pbcn, (d) C2/c, (e) Cmca-12, and (f) Cmca-A at 
250 GPa. The right hand panels show the vibrational densities of states. The modes shown with negative frequencies indicate 
dynamical instabilities. The modes below 2500 cm -1 are lattice modes and those at higher frequencies are vibronic modes. 
The densities of states of the Pc, Pc-96 and Pbcn mixed phases are very similar and each has three distinct bands of vibrons. 
Pc has some unstable modes near Z while the unstable modes of Pbcn cover a substantial fraction of the zone. The Pc-96, 
C2/c, Cmca-12 and Cmca-A structures are dynamically stable. 
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FIG. 4. (color online). Enthalpies/free energies of the Cmca- 
12, Pc and Cmca-i structures relative to C2/c for a) the static 
lattice structures, b) with proton ZP motion and c) with full 
vibrational motion at 300 K. The static-lattice enthalpy of 
Cmca-4: is too high for it to appear in a) . 
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FIG. 5. (color online). Computed phase diagram of hydrogen 
at high pressures. 



pressure. The pressure dependence of the Raman vi- 
bronic frequency of C2/c matches that observed for phase 
III reasonably well. The Cmca-4 structure^ is not con- 
sistent with phase IV as it has a single intense Raman 
vibron mode at frequencies of about 3000 cm -1 which 
slowly declines with pressure and, as mentioned above, 
Cmca-4 is significantly less stable than other phases at 
250 GPa. Pc has two strong vibronic Raman modes, one 
at just under 4200 cm -1 whose frequency is very weakly 
pressure dependent, and a second and more intense vi- 
bronic mode at lower frequencies which softens rapidly 
with increasing pressure, see Figs. [6] and [7] and Table fl] 
Inspection of the vibrational eigenmodes shows that the 
higher frequency mode arises from the strongly bonded 
molecular layers while the lower frequency mode arises 
from the weakly bonded layers. Eremets and TroyaiJ^ 
did not discuss the less-intense higher-frequency Raman 
vibron mode of phase IV, but it is visible in Fig. 2a of 
their paper and in Fig. S3 of their supplementary infor- 



mation, and in the Raman data of Figs. 2 and 3 of Rcf. 
IT5I . The higher frequency vibron of the Pc structure is a 
good match for the Raman mode of phase IV observed at 
about 4130 cm" 1 in the range ~220-250 GPa (see Fig. 
2a of Ref. [lil ). The higher experimental Raman vibron 
frequency and that calculated for Pc are almost inde- 
pendent of pressure in this range, see Fig. [S] and Table 
III 

The experimental room temperature Raman data for 
hydrogen in Fig. 3 of Ref. [ID shows phase IV forming 
from phase III as the pressure is increased, see Fig. [71 
which gives rise to the "knee" in the lower frequency vi- 
bron branch at around 220 GPa. The downward slope 
of the room temperature experimental data above 220 
GP a 14 i 15 is in reasonable agreement with the calculated 
slope of the Raman data for the lower vibron mode of Pc, 
although the calculated modes are about 500 cm" 1 lower 
in frequency, see Figs. [5] and [7] and Table Q] The neglect 
of the effects of proton motion on the vibrational modes 
in our calculations might have a significant effect on the 
frequencies, as might the approximate density functional. 
The Pc structure does, however, yield qualitative agree- 
ment with the two Raman vibrons observed at room tem- 
perature in phase I V 14 i 15 , whereas the C2/c, Cmca-12 
and Cmca-4 structures provide very poor descriptions. 

The lattice Raman modes of C2/c, Cmca-12 and Pc 
extend to frequencies just above 2000 cm -1 at 250 GPa. 
Pc shows substantially higher lattice-mode Raman in- 
tensity than C2/c or Cmca-12. C2/c does not have any 
very strong lattice Raman modes, but Cmca-12 has a 
fairly intense mode at around 1100 cm -1 . The strongest 
peak in the lattice Raman modes of Pc is at 360 cm -1 at 
250 GPa, with other peaks at 110, 730, and 1470 cm" 1 . 
Phase IV shows three strong lattice-mode peaks at 250 
GPai£, although peaks below about 150 cm -1 might be 
obscured by the filter. The strongest peak observed at 
300 enr 1 is in reasonable agreement with our peak at 360 
cm , and the other two experimentally observed modes 
at 630 and 1065 cm -1 are in rough agreement with our 
calculated values of 730 and 1470 cm -1 . We find qual- 
itative agreement between our calculated lattice Raman 
modes for Pc and those observed in phase IV—, whereas 
the C2/c, Cmca-12 and Cmca-4 structures provide very 
poor descriptions. 



F. Infra-red spectra 

We have also calculated IR spectra for C2/c, Cmca-12 
and Pc, see Fig. [5] and Table Q] The vibronic modes of 
C2/c and Cmca-12 soften with increasing pressure. C2/c 
has a single intense IR vibron whose frequency drops with 
pressure from 4150 cm" 1 at 200 GPa to 4030 cm" 1 at 300 
GPa. Cmca-12 has two intense IR vibrons at somewhat 
lower frequencies than the single intense vibronic mode 
of C2/c. The frequencies of the two intense IR vibrons 
of Cmca-12 arc about 3790 and 4030 cm" 1 at 200 GPa, 
and these frequencies fall to 3670 and 3940 cm -1 , re- 
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FIG. 6. (color online). Raman spectra of C2/c, Cmca-12 and 
Pc at a) 200 GPa and b) 250 GPa. 
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and Pc as pressure increases. The IR vibron of C2/c is 
considerably more intense than those of Pc, and therefore 
we expect a drop in IR vibron intensity at the transition 
from phase III to IV. 

As shown in Fig.[5]and TablelH the IR lattice modes of 
C2/c, Cmca-12 and Pc stiffen with increasing pressure. 
C2/c has a single intense lattice mode whose frequency 
of 1830 cm" 1 at 200 GPa increases to 2120 cm" 1 at 300 
GPa. This mode is in reasonable correspondence with 
observations for phase III of hydrogen. Cmca-12 has 
three intense lattice modes at 1380, 1730 and 1940 cm -1 
at 200 GPa, which stiffen to 1550, 1960 and 2230 cm" 1 
at 300 GPa. Pc shows two intense lattice modes at 1150 
and 1850 cm -1 which stiffen to 1250 and 2110 cm -1 at 
300 GPa. 

IR data for hydrogen up to 360 GPa and 250 K has 
recently been published—. According to our phase dia- 
gram (Fig. [5]) none of the data sets in the lower part of 
Fig. 3 of Ref. [H are within the region of stability of Pc, 
but the lower pressure data are within the region of C2 / c 
and the higher pressures data at (330 GPa, 100 K) and 
(360 GPa, 17 K) are within the region of Cmca-12. The 
drop in the observed IR vibron frequency visible in Fig. 
3 of Ref. [H of about 240 cm -1 between the data at (292 
GPa, 200 K) and (330 GPa, 100 K) is consistent with a 
phase transition from C2/c to Cmca-12, if the broaden- 
ing of the peaks at 330 GPa and 360 GPa is taken to arise 
from the unresolved pair of vibronic peaks of Cmca-12, 
whose mean frequency is about 230 cm -1 lower than the 
vibron frequency of C2/c at 300 GPa. 
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FIG. 7. (color online). Pressure dependence of the Raman 
vibron frequencies of C2/c (blue solid line), Cmca-12 (red 
solid line) and Pc (black solid line). The experimental data 
for phase III at 100 K (dashed orange line) are from Ref. 03 
and the room temperature data for phases III and IV are from 
Ref. [l5l (dashed green line). 



FIG. 8. (color online) Infra-Red spectra of C2/c, Cmca-12 
and Pc at a) 200, b) 250, and c) 300 GPa. 



G. Bandstructures 



spectively, at 300 GPa. Pc has two intense IR vibrons 
at 3760 and 4520 cm" 1 at 200 GPa, but the intensity of 
the lower frequency mode falls rapidly with pressure, and 
at 300 GPa Pc has only a single intense IR vibron peak 
at 4600 cm -1 . The higher frequency IR vibron of Pc 
stiffens with pressure. The strength of the Cmca-12 IR 
vibron modes increases in comparison to those of C2/c 



The DFT bandstructures of C2/c, Cmca-12 and Pc at 
250 GPa are shown in Fig. [9j and the pressure depen- 
dence of the band gaps is shown in Fig. [101 Cmca-12 
has the smallest gap and Pc has the largest. The ther- 
mal (minimum) band gaps decrease with pressure and 
go to zero at 210 GPa (Cmca-12), 260 GPa (C2/c), and 
320 GPa (Pc), while the optical band gaps (minimum 
direct gaps) are substantially larger. Standard DFT ap- 
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TABLE I. Frequencies in cm 1 of the main IR and Raman peaks of C2/c, Cmca-Yl and Pc as a function of pressure. 
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proaches, however, normally underestimate band gaps, 
and exact-exchange DFT calculations indicate that the 
true gaps may be roughly 1 eV larger—. Increasing the 
band gaps by 1 eV gives metallization pressures of about 
250 GPa (Cmca-12), 310 GPa (C2/c), and 370 GPa (Pc). 
We have neglected the effect of nuclear vibrational mo- 
tion on the band gaps, which is expected to be significant. 
The inclusion of vibrational motion is likely to reduce the 
band gap because it amounts to averaging the gap over 
the structures accessed during the motion, which are less 
stable than the static lattice structure and therefore tend 
to have smaller band gaps. This effect is expected to be 
significant in hydrogen because of the large vibrational 
amplitudes, and it may be particularly important in the 
weakly bonded layers of the Pc phase. The underestima- 
tion of the band gaps arising from the PBE functional 
may tend to be cancelled by our neglect of the effects 
of coupling between the electronic and vibrational de- 
grees of freedom. We cannot therefore obtain a reliable 
estimate of the gap at this level of theory but, if we sup- 
pose that the DFT calculations underestimate the band 
gap by 1 eV, we find that the coupling of the electronic 
and vibrational motion would have to reduce the gap by 
about 1 eV to reproduce the room-temperature gap of 
phase IV of 1.8 eV at 315 GPa reported by Howie et 
al^. Theoretical work describing the anharmonic pro- 
ton vibrations at finite temperatures and their effect on 
the structure and band gap will be required to investigate 
this in more detail. 



IV. SUMMARY AND CONCLUSIONS 

In summary, our results support the observation by 
Eremets and Troyanii and Howie et al^- of a structural 



phase IV of hydrogen at room temperature and pressures 
above 200 GPa. The C2/c, Cmca-12 and Cmca-A candi- 
date structures cannot account for the observed proper- 
ties of phase IV. DFT calculations predict the Pc struc- 
ture to be stable at room temperature over the pressure 
range 250-295 GPa, in reasonable agreement with the 
experimental observations of phase IV. Mixed structures 
are stabilized by the presence of low-frequency vibra- 
tional modes arising from the weakly-bonded graphenc- 
like layers. Mixed structures explain the observation of 
both the high frequency Raman vibron of phase IV that 
is almost independent of pressure and arises from the 
short bonds in the strongly-bonded layers, and the lower 
frequency vibron whose frequency falls rapidly with pres- 
sure, which arises from the weakly bonded layers. We 
find the C2/c structure, that was previously identified 
as a candidate for phase IIli£, to be stable at low tem- 
peratures. The larger lattice-mode Raman intensity of 
Pc compared with that of C2/c is also consistent with 
experiment, and we find qualitative agreement between 
the calculated lattice Raman modes of Pc and those ob- 
served for phase I V 14 i 15 . We conclude that mixed phases 
are reasonable structural models of phase IV of hydrogen. 
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